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We review the recent developments in the field of photonic lattices emphasizing their 
unique properties for controlling linear and nonlinear propagation of light. We draw some 
important links between optical lattices and photonic crystals pointing towards practical 
applications in optical communications and computing, beam shaping, and bio-sensing. 
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1. Optical lattices as nonlinear photonic crystals 

Nonlinear propagation of light in periodic structures has become an attractive area 
of research in recent yearsSMH, holding strong promises for novel applications in 
photonics. The underlying physical effects are analogous to those occurring in a 
number of different systems, including biological molecular structures, solid-state 
systems, and Bose-Einstein condensate in periodic potentials. A special class of pe- 
riodic structures are photonic crystals 4 ", artificial materials with optical bandgaps, 
which offer unique possibilities for controlling the propagation of light in a way 
similar to that of semiconductors used for manipulating the flow of electrons. Pho- 
tonic crystals represent a broad class of structures with periodicity of the refractive 
index on the wavelength scale in one, two, or three dimensions. They were first 
suggested in the pioneering papers of Eli Yablonovitch 5 and Sajeev JohiP, and 
nowadays this terminology is widely applicable to many different materials (see 
Fig. [T]) , some of which are commonly used in almost any optical laboratory. Exam- 
ples of one-dimensional photonic crystals include dielectric mirrors, Bragg gratings, 
and arrays of optical waveguides3H Two-dimensional photonic crystals are com- 
monly represented by photonic crystal fibers^ and planar photonic crystals^. In 
three dimensions, photonic crystals are known to exist in opal, wood-pile, or in- 
verse opal geometries^. The common feature of all these different structures is that 
they allow for manipulation of the flow of light in the direction of periodicity. There- 
fore, photonic crystals offer a possibility to achieve ultimate control over linear and 
nonlinear properties of the light propagation, as well as enhanced control over light 
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Fig. 1. Examples of one-, two-, and three-dimensional photonic crystals 



emission and amplification. Herewith we concentrate primarily on the abilities of 
the periodic structures to control both linear and nonlinear propagation of light. 

Several approaches for controlling the light propagation by engineered periodic 
structures have been suggested theoretically and demonstrated experimentally in 
recent years. These include manipulation of linear light propagation (refraction, 
diffraction, and dispersion) as well as various nonlinear effects such as harmonic 
generation, stimulated scattering, and nonlinear self-action. Important examples 
for a control of refraction and diffraction of light constitute the effects of negative 
refraction^ and self-collimation 11 . Furthermore, a design of the structural disper- 
sion allows one to manipulate the group velocity of light resulting in the slow-light 
propagation^, or enhancement of the spectral response with a fascinating example 
of the superprism effect 13 . A nonlinear response of the material offers novel oppor- 
tunities for dynamic tunability of the structures by varying the light intensit}^. An 
interplay between nonlinearity and periodicity represents unique way to efficiently 
manipulate light by light for optical switching and signal processing applications. 
Furthermore, light control can be scaled down to micron-scale structures suitable 
for integration of multiple functionality on a photonic chip. 

The goal of our studies in this field is twofold. First, we analyze the fundamentals 
of nonlinear wave physics in periodic media and, second, through relevant exper- 
imental observations we exploit these effects for practical applications. Currently, 
the major research on photonic structures is concentrated on the demonstration of 
fundamental physical phenomena. However, novel photonic applications are show- 
ing up on the horizon promising a significant impact on various technologies in 
optical communications and computing, beam steering and shaping, bio-sensing 
and medical diagnostics. 

To achieve our goals we need to employ periodic structures of different geome- 
tries, which are easy to fabricate, and which possess strong nonlinearities at mod- 
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erate laser powers. Currently there exists a number of different approaches for the 
fabrication of photonic crystals including focused ion beam milling, e-beam lithog- 
raphy combined with reactive ion etching, two-photon and UV polymerization. All 
these techniques, however, are resource demanding and cost ineffective, imposing 
great constrains on the fundamental research, where flexible and quick modification 
of the structural parameters is required. A simpler fabrication and characterization 
of the periodic photonic structures can be achieved when the scale of periodic- 
ity is larger than the wavelength of light. Termed as optical lattices such periodic 
photonic structures include optical waveguide arrays^, photonic crystal fibers^, and 
optically-induced lattice3l l l 17 l 18 [ The scale of periodicity in such structures is of the 
order of a few micrometers, and their fabrication is facilitated by standard litho- 
graphic techniques, fiber drawing, or multiple beam interference. The challenge is 
to achieve strong nonlinear response of the material at moderate laser powers. The 
nonlinear response can be enhanced through extended propagation length (as in the 
case of fibers), stronger light confinement (as in nanowires and photonic crystals), 
or through slow and resonant nonlinearities (as in the case of optically-induced 
lattices). In our current studies we chose materials with slow nonlinear response, 
but with access to nonlinear effects at micro- Watt laser powers. The periodic pho- 
tonic structures with strong nonlinear response and various geometries allow us to 
study the basics of nonlinear physics in periodic structures. The obtained knowledge 
can then be applied to the fabricated photonic structures and photonic crystals in 
particular, thus bringing the advantages of miniaturization and integration on a 
photonic chip. We believe that the process of basic knowledge accumulation and its 
further application is the successful path towards future photonic technologies. 

In this paper, we review some recent advances and fundamental concepts of 
nonlinear light propagation in periodic photonic structures, emphasizing their abil- 
ities to control the spatial dynamics of light propagation. The paper is organized as 
follows. In Section |2j we present the basic concepts of the optically-induced lattices 
and discuss how their tunability can be used to engineer the linear light propaga- 
tion such as refraction and diffraction. In Section |3j we review our results on the 
nonlinear localization of light in periodic structures and outline some novel ideas for 
the nonlinear control of wave transport in periodic structures. Section [4] is focused 
on novel physical phenomena of light propagation in fabricated periodic photonic 
structures, including the beam interactions with the interfaces. 

2. Fundamentals of light propagation in periodic structures 

The propagation of light in periodic optical lattices is determined by three major 
effects: (i) the coupling between neighboring sites of the lattice (inter-site coupling), 
(ii) Bragg scattering arising from the periodicity of the lattice, and (iii) the spe- 
cific lattice geometry, being one- (ID), two- (2D), or three-dimensional (3D). The 
geometry is of particular importance in 2D and 3D, where the lattice symmetry 
(e.g. square, hexagonal) significantly affects both the inter-site coupling and wave 
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scattering in different directions. Being able to engineer all those characteristics, 
one can fully control the propagation of light. Additionally, dynamic tunability 
of the lattice depth and periodicity in real time will provide an ideal system for 
experimental studies of fundamental physical phenomena. 



2.1. Optically-induced one- and two-dimensional lattices 

A great opportunity for controlling the lattice parameters was offered by the idea of 
optically-induced lattices in a biased photorefractive crystal proposed by Efremidis 
and co-workers 16 . With this theoretical proposal, the optical lattice is induced by 
the interference of two or more broad laser beams propagating inside the photore- 
fractive crystal. The authors took advantage of the strong electro-optic anisotropy 
of the crystal, such that the ordinarily polarized lattice-forming beams are not af- 
fected by the applied electric field. The strong but anisotropic electro-optic effect of 
the crystal leads to index changes predominantly for the extraordinarily polarized 
light. Thus, the periodic light pattern resulting from the multiple beam interfer- 
ence will induce a periodic optical potential for any extraordinarily polarized probe 
beam. At the same time, probe beams will experience strong nonlinear self-action 
at moderate (/iW) laser powers. The propagation of light in such a system can be 
described (in isotropic approximation) by the nonlinear Schrodinger equation for 
the slowly varying amplitude of the electric fiel d 16 I 19 I22I 

dE ^ (d 2 E d 2 E\ 7 V / \ 

where (x,y) and z are the transverse and propagation coordinates normalized 
to the characteristic values x s = y s = 1 fim and z s = 1 mm, respectively, 
D = z s \/ (Aim^x 2 ) is the diffraction coefficient, A is the wavelength in vacuum, 
no is the average refractive index of the medium, and Vq is the bias voltage applied 
across the crystal. The term jVo(Ib-\-I p (x, y) + \E\ 2 )~ 1 characterizes the total refrac- 
tive index modulation induced by the optical lattice and the probe beam, linearly 
proportional to the applied electric field. Here = 1 is the normalized constant 
dark irradiance, and the lattice intensity I p depends on the specific lattice geometry 
controlled by the number and position of the lattice forming beams. For example, 
I p (x,y) = I g \ex.p(ikx) + exp(— ikx)\ 2 for ID lattice produced by two-beam inter- 
ference [Fig. [5Ja)]; I p (x,y) = I g \ ex.p(ikx) + exp(— ikx) + ex.p(iky) + exp(— iky)\ 2 
for a 2D lattice of square geometry created by four coherent beams [Fig. [2^b)] ; 
or I p (x,y) = I g \exp(ikx) + exp(— ikx/ 2 — iky\fZ/2) + exp(— ikx/ 2 + iky V3/2)\ 2 
in the case of a hexagonal lattice generated by three-wave beam interference 
[Fig. fflc)]. Other configurations of optical patterns can be used to induce 
BesseP 1 ! 2 ^^ 3 !^ 4 !^ 5 !^^!^ 7 !^ 8 !^ 9 ! 30 !!!!!!!, Mathieu 33 , and quasi-periodic 34 lattices. 

The advantage of using optically induced lattices as a test-bed physical sys- 
tem comes from the fact that the inter-site coupling can be easily controlled via 
the applied bias voltage Vq (positive or negative) or lattice-forming beam intensity, 
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Fig. 2. Optical induction of (a) one- and (b,c) two-dimensional photonic lattice in a biased 
photorefractive SBN crystal of square and hexagonal lattice, respectively. 

while the Bragg scattering can be controlled via the period of the lattice, which 
is inversely proportional to the angle between the interfering beams. A disadvan- 
tage of the system is the relatively low index modulation, typically of the order of 
few times 10 -4 . This disadvantage results in the need of a rather large lattice pe- 
riod, and relatively short propagation distances (typical crystal is 1-2 cm long). The 
possibility for dynamic reconfiguration of the lattice, however, makes the optically 
induced lattices an attractive tool for studies of nonlinear phenomena in periodic 
photonic structures of various geometries. Indeed, this technique was applied by 
several groups and the experimental demonstration of various fundamental effects 
came 



2.2. Bandgap structure and tunability 

Even though the technique of optical induction results in weak index modulation, 
this can be sufficient for the appearance of a distinct bandgap structure for the 
transverse components of the wave vectors K x1 K y due to the Bragg scattering 
of waves propagating at small angles. The propagation of linear waves through a 
periodic lattice is described by the spatially extended eigenmodes, called Bloch 



waves. They can be found as solutions of linearized Eq. (2.1) in the form E = 
ip(x, y) ex.-p(i/3z + iK x x + iK y y), where ip(x, y) has the periodicity of the underlying 
lattice, and (3 is the propagation constant. 

When the applied bias voltage is zero, the photorefractive crystal is homo- 
geneous and the dispersion relation for the propagation of waves is parabolic, 
(3 = —D(K^ + Ky), as shown in Fig. [3|a). Once a voltage is applied across the 
crystal, the dispersion is modified, and there appear forbidden gaps in the trans- 
mission spectrurri 7 ! 37 !, as shown in [Fig. [3^b)] for the case of ID lattice. Note that 
the dispersion curves become periodic, /3(K X ) = /3(K X + 27r/d), where d is the 
lattice period, and accordingly are fully characterized by their values in the first 
Brillouin zone, —n/d < K x <ir/d. The ID band-gaps open around the points with 
K x d/7T — 0, ±1, ±2, . . ., where the Bragg scattering condition is satisfied. 

The modification of the dispersion relation dramatically changes the beam prop- 
agation. The direction of propagation is determined by the normal to the dispersion 
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Fig. 3. (a) Dispersion of linear waves in a bulk medium, and (b,c) Bloch-wave dispersion in 

(b) weaker and (c) stronger ID optical lattices. Arrows indicate propagation directions at selected 
dispersion points K x = 7r/(2d). 



curves [see the arrows in Figs. [3]^b),(c)], while diffraction is determined by the cur- 
vature at the corresponding point. Thus waves associated with the top of the first 
band and propagating along the lattice will experience normal diffraction, while 
waves corresponding to the bottom of the first band (convex curvature) will expe- 
rience anomalous diffraction 38 . The curvature of the dispersion curves changes in 
between those two zones from concave to convex, therefore there will be a point 
where the curvature is zero. Waves belonging to this zone will, therefore, experi- 
ence no second-order diffraction, similar to the self-collimation effect in photonic 
crystaleP. 

It is important to note that the size of the forbidden gaps is proportional to 
the induced index modulation. An increase of the index modulation leads to flat- 
tening of the bands and therefore a wider Bragg reflection gap. Therefore, if one 
increases the index modulation by increasing the applied bias field, one can dy- 
namically control the direction of propagation of waves inside the structure. As 
illustrated in Figs. J3^b) and (c), the increase of the applied voltages makes the 
waves with a normalized transverse momentum of 0.25 to propagate under smaller 
and smaller angles inside the structure. Such dynamic tuning of wave propagation 
is not an isolated phenomenon for optically-induced lattices, but can be applied 
to other types of periodic structures with electro-optic tunability. In particular, 
various possibilities for beam control were demonstrated in liquid-crystal cells with 
periodical modulation induced by an array of electrode j 4Q | 41 | 42 | 41 | 43 | 44 | 4 §l 

The optically induced bandgap structure is also well defined in the case of 
periodic lattices with higher dimensionality. In 2D lattices, the band-gaps may 
appear only above the critical value of the lattice depth. Examples of characteristic 
bandgap structures for square and triangular lattices are presented in Fig. [4] Note 
that a complete 2D bandgap appears between the first and the second band. This 
fact underlines the importance of the optically-induced lattices as an analogue of 
ID and 2D photonic crystals. Therefore, the optically induced lattices represent an 
accessible test-bed for studies of generic bandgap phenomena in photonic periodic 
structures. 
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Fig. 4. Bloch-wave dispersion in 2D optical lattices of (a) square and (b) triangular symmetries. 
The horizontal axes correspond to the contours passing through the high-symmetry points of the 
Brillouin zones shown in the inset. 

2.3. Selective excitation of Block waves 

In order to explore the ability for control and steering of beams in an optical lat- 
tice, it is important to understand the character and the profiles of the propagating 
linear waves. Finite beams propagating inside the lattice can be naturally repre- 
sented as a superposition of Bloch waves of the induced periodic potential. The 
Bloch wave profiles corresponding to the top and bottom of the first band and the 
top of the second band of the ID lattice bandgap spectrum [Fig. [5|a)] are given 
in Fig. [5|b, solid line), together with the underlying periodic index modulation 
Fig. [5jlo, shading). These Bloch waves have the periodicity of the lattice itself, but 
can also have a nontrivial phase structure. The Bloch wave from the top of the 
first band is periodic, but always sign-definite with a fiat phase front. On the other 
hand, the Bloch waves from the bottom of the first band and the top of the second 
band change sign at each lattice period. Furthermore, the maxima of the Bloch 
wave intensity associated with the bottom of the first band are centered on the 
maxima of the refractive index modulation, whereas the intensity maxima of the 
Bloch waves corresponding to the top of the second band are centered in between 
the lattice maxima. 

Each Bloch wave has a well defined propagation direction and diffraction coef- 
ficient, following the dispersion curves of the bandgap diagram. In order to explore 
their specific properties, it is important to know how to selectively excite each Bloch 
wave. Experimentally, this means that most of the light entering the lattice should 
be efficiently coupled to the desired Bloch wave. For best excitation one needs to 
exactly match the field profile of the corresponding Bloch wave. Exact matching of 
the Bloch wave profile, however, is not always possible. A good approximation is 
given by the leading-order Fourier component of the Bloch wave profile, plotted in 
Fig. |5|b) with a dashed line. The leading-order Fourier component for the Bloch 
wave associated with the top of the first band is a constant beam intensity, while 
for the Bloch waves from the edges of the Bragg reflection gap, the leading order 
Fourier component is cosine or sine functions. Having this in mind, we demonstrated 
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Fig. 5. Bloch waves and their excitation: (a) typical bandgap diagram for ID optical lattice, (b) 
Profiles of the Bloch waves (solid line) corresponding to the edges of the first and second bands, 
superimposed on the leading-order Fourier component (dashed line) and induced refractive index 
change (shading). (c,d) Experimental excitation of Bloch waves: input and output of the lattice, 
respectively. Insets show the excitation geometry. Dashed line - output beam profile without the 
lattice. 



experimentally efficient excitation of the Bloch waves by a broad Gaussian beam 
for the Bloch wave of the top of the first band and by two-beam interference for 
the Bloch waves from the bottom of the first band and top of the second band, 
respective!}®. The only difference in the excitation between the latter two is the 
position of the interference pattern with respect to the lattice [see Fig.[5jc)], being 
centered at or in between the lattice sites. The position of the interference maxima is 
easily controlled by the relative phase between the two input beams (see the bottom 
inset in Fig. [5|. The experimental profiles of the excited Bloch waves at the crystal 
output are shown in Fig. |5|d) and clearly demonstrate that our technique leads to 
efficient excitation of the corresponding Bloch waves. This is particularly evident 
in the excitation of the Bloch waves from the top of the second band which have 
a characteristic double peak structure, successfully reproduced in our experiments. 
Other approaches for selective Bloch- wave excitation were also demonstrated based 
on wave- vector matching with side-on beam incidence^ or prism coupling^. Note 
that if the input is not matched to a particular mode, then a spectrum of Bloch 
waves can be excited simultaneously from several bands, demonstrating complex 
interactions in the nonlinear regimd 49 l 50 l 51 l 52 l 53 l 54 l 55 l 56 l 57 l 58 T^H 

We have also realized selective Bloch- wave excitation in 2D lattices of squardHEH 
and triangular^ symmetries. The calculated Bloch- wave intensity and phase struc- 
ture for the high symmetry points of the 2D square lattice from the first and 
second spectral bands are shown in Fig. [6] The upper row shows the Bloch-wave 
intensity profiles and the next row shows the corresponding phase structure. As 
a reference, the first column shows the light intensity of the lattice itself and the 
corresponding Brillouin zone. For the two-dimensional Bloch waves from the first 
band, the intensity distribution of all modes reflects the structure of the square 
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Fig. 6. 2D Bloch waves and their excitation: Intensity (top row) and phase (second row) of 
different Bloch modes from the high symmetry points of the first and second bands of a square 
lattice. The blue color for the phase distribution corresponds to the zero phase, while the red 
color corresponds to the 7r phase. Bottom rows: experimental intensity profiles for excitation of 
2D Bloch modes at the input (third row) and the output (bottom row) of the lattice. 



lattice, with the intensity maxima coinciding with those of the lattice. However, 
their phase structures differ substantially. As can be seen from Fig. [6j the phase 
of the two-dimensional Bloch waves originating from the I\ point is constant. The 
phase structure becomes nontrivial for the modes from the Xi and Mi points. For 
the Xi (Yi) point, the phase represents a stripe-like pattern being constant along 
one principal direction of the lattice and exhibiting tt phase jumps along the other 
direction. For the Bloch waves originating from the Mi point, the phase distribu- 
tion resembles a chessboard pattern. On the other hand, the two-dimensional Bloch 
modes from the second spectral band have the intensity maxima centered between 
the maxima of the square lattice. The phase structure has a form of stripes oriented 
along one of the principal directions of the two-dimensional lattice for the X2 point, 
or in 45° with respect to the principal axes for the M2 point. 

The difference in the phase structure of the two-dimensional Bloch waves trans- 
lates into differences in the propagation dynamics for beams of a finite size which 
spectrum is localized in the vicinity of the corresponding high-symmetry points in 
the Brillouin zone. Indeed, the alternating phase is a signature of strong Bragg 
scattering, that may lead to an enhanced diffraction of beams similar to the ef- 
fect of the dispersion enhancement in the Bragg gratings^. Therefore, the beams 
can experience anisotropic diffraction due to the asymmetric phase structure of 
the corresponding Bloch waves, and this can be detected by analyzing the beam 
broadening in the linear regime. 
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Fig. 7. Tunable refraction in a tilted lattice: (a) Scheme of beam propagation, (b) Tuning of 
beam refraction for the different Bloch modes by varying the applied bias field. The lattice tilt 
is 2 mrad. (c) Tunable refraction by varying the lattice tilt at a constant bias field 3kV/cm. (d) 
Corresponding output intensity profiles for a lattice tilt of 3 mrad. 



In order to study experimentally the generation, formation, and propagation of 
Bloch waves in two-dimensional photonic lattices, we have developed a special setup 
to match the structure of the input optical beam and the selected Bloch wave, see 
the third row in Fig. [6] This is achieved by the use of the programmable phase mod- 
ulator that converts the initially Gaussian probe beam into the desired amplitude 
and phase modulation at the front face of the photorefractive crystal. A finite beam 
will diffracts with a rate depending on the value of the diffraction coefficients along 
the principal directions of the lattice. This allows one to characterize the lattice 
dispersion by analyzing the diffracted beam patterns shown in Fig. [6j bottom row. 



2.4. Tunable refraction 

The efficient excitation of Bloch waves can be used in experimental schemes for 
steering and switching of beams in periodic structures. For example, coupling of 
light into a particular Bloch wave can dramatically change the direction of beam 
propagation and thus the transport of energy inside the structure can be controlled. 
To demonstrate this idea, we performed experiments on tunable beam refraction 
and beam steering in optically induced lattices^. In a straight lattice, all Bloch 
waves from the top and bottom of each band propagate exactly along the lattice 
(zero transverse wavevector component). The different Bloch waves, however, can 
be separated if the lattice is tilted at a small angle as shown in Fig. (7^a). Such a 
lattice tilt translates into a tilt of the bangap structure, causing Bloch waves from 
the first band to bend (or refract) along the tilt of the lattice, while Bloch waves 
from the second band bend in the opposite direction [Fig. [T^a)]. The change in the 
direction of propagation leads to an output beam shift if the initial light is coupled 
to a specific Bloch wave. 

As already discussed, the direction of propagation of each Bloch wave will de- 
pend on the depth of the induced refractive index modulation and therefore the 
output shift can be also tuned by changing the bias voltage applied to the crystal. 
Experimentally, we measured this shift as a function of the applied voltage across 
our 5 mm thick crystal. Our experimental results are depicted in Fig.[7Fb). It is clear 
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from the figure that with increasing index modulation (bias voltage ~ 2.5 kV), all 
the Bloch waves approach the tilt of the lattice given by the dashed line in the 
figure. At lower voltages (1 — 2kV), however, we observe that Bloch waves from the 
top of the first band are shifted less than the lattice, while those from the bottom 
of the first band and top of the second band are shifted substantially more. We 
measured a six-fold increase of the lattice shift (at lkV) with positive and neg- 
ative gain for Bloch waves associated with the bottom of the first band and top 
of the second band, respectively. These experimental observations demonstrate an 
important example of tunable refraction in periodic structures, pointing towards 
potential applications in beam steering technologies. 

Furthermore, we measured the dependence of the beam shift at the output 
vs. the lattice tilt [Fig. [7jc)]. Our results showed that for a small lattice tilt, the 
dependence is linear, while it saturates to a value equal to the Bragg angle inside 
the lattice for Bloch waves from the bottom of the first band and top of the second 
band. The corresponding output beam profiles of the three different Bloch waves 
are shown in Fig. [7^d). They show that in each case the input light is coupled 
directly to a particular Bloch wave, with no mixed excitation, and that the beams 
experience a strong transverse shift. 

3. Nonlinear localization and gap solitons 

In order to confine optical beams inside the periodic structure it is necessary to 
balance the beam diffraction through nonlinear self-action. However, as we already 
know, the diffraction of beams can be significantly modified by the effects of period- 
icity. Beams associated with the top of each band have normal diffraction which can 
be balanced by a self-focusing type nonlinearity. Beams from the bottom of each 
band, on the other hand, exhibit anomalous diffraction, which can be balanced by 
a self-defocusing type nonlinearity. 

3.1. Solitons in periodic structures 

In a bulk medium with Kerr-type nonlinearity, the propagation of light of high 
intensity will result in a local refractive index change of the material. For focus- 
ing nonlinearity, this index change is equivalent to an induced optical waveguide 
[Fig. [8^a, top)] which can trap and guide the beam in the structure, resulting in 
the formation of spatial solitons^ [Fig. |8|b, top)]. In the case of defocusing nonlin- 
earity, the induced index change is negative [Fig. [8^a, bottom)] and will result in 
anti-guiding of the light, or beam self-defocusing. Therefore, no localized ("bright") 
solitons can exist in the case of bulk defocusing material [Fig. [8^b, bottom)]. This 
picture, however, changes fundamentally in periodic structures. In the case of focus- 
ing nonlinearity solitons exist near the top of each band, where diffraction is normal. 
This can result in the formation of discrete spatial solitons from the top of the first 
band 66 . Discrete solitons can be excited by a narrow input beam launched into a 
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Fig. 8. Formation of spatial solitons: (a) self-focusing and self-defocusing nonlinearity will induce 
a positive or negative defect in the material, respectively, (b) In a bulk material solitons can form 
only for self-focusing nonlinearity. In a periodic medium solitons can form from the first (c) and 
second bands (d) for both focusing and defocusing nonlinearity. 



single lattice site [Fig. [8jc, top)]. They were first observed in AlGaAs waveguide 
arrays^HHl an( j } a ^ er reproduced in optically induced latticedMlS, and liquid crystal 
waveguide arrays^. In the case of higher order bands, spatial solitons have been 
predicted to exist in the Bragg reflection gaps' 70 ! 71 ! Their excitation is less trivial, 
and both side-on excitation^ and head-on excitation with periodically modulated 
input beamdHES have been demonstrated [Fig. [8^d, top)]. The latter appears ben- 
eficial when the transverse soliton velocity needs to be controlled, including the 
special case of excitation of immobile gap solitons. 

In the case of defocusing nonlinearity ID solitons exist near the bottom of each 
band and their propagation constant moves down into the Bragg reflection gap 
with increasing beam intensity. Solitons in periodic waveguide arrays with defocus- 
ing nonlinearity were already predicted in 1993^. They possess a staggered phase 
structure and the propagation constant lies inside the first Bragg reflection gap 
[Fig. |8](c, bottom)]. Therefore, such solitons have properties similar to both dis- 
crete and gap solitons in self- focusing nonlinear materials. Herewith we will refer 
to them as staggered solitons. Staggered solitons were first demonstrated experi- 
mentally by excitation with an inclined input beanJUEH With this method, only 
a fraction of the input light is coupled into the staggered soliton mode. Later, 
the generation of staggered solitons was demonstrated in defocusing lithium nio- 
bate waveguide arrays by excitation with a TEMio laser moddSEIl Such excitation 
leads to coupling to the unstable "even" soliton mode, which then transforms to the 
stable "odd" -symmetry mode. The similarity of the staggered solitons and discrete 
solitons, however, points to a simpler excitation method of high efficiency. As we 
recently demonstrated 78 efficient excitation of staggered solitons can be achieved 
with a narrow input beam coupled to a single lattice site [see Fig. [8^c, bottom)], 
provided the refractive index contrast of the structure is high enough. In this case 
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the beam propagation inside the lattice is very similar to the propagation in lat- 
tices with focusing nonlinearity, with the only difference being the staggered phase 
structure of the output beam. Therefore, one can see that if the optical lattice can 
be treated as a discrete system, there is a direct analogy between the localization 
of light for focusing and defocusing nonlinear it ies. 

3.2. Two-dimensional solitons and their mobility 

The localization of optical beams is even more interesting and nontrivial in the case 
of higher dimensions. In a 2D lattice, for example, waves can experience anisotropic 
diffraction 7 - 9 -^, which strongly affects the formation of solitons. For a square lattice 
with a period of 23 /am [Fig.[9^a)] optically induced in a photorefractive SBN crystal, 
the typical bandgap structure folded along a contour through the high-symmetry 
points of the lattice is shown in Fig. [4^a). For self- focusing nonlinearity, solitons 
are associated with the top of each band, where diffraction is normal. Due to the 
crossing of higher order bands, there are only two complete gaps in our 2D lattice: 
the total internal reflection gap and one Bragg reflection gap. This is in sharp con- 
trast to the ID case, where multiple Bragg reflection gaps usually exist. Therefore, 
the soliton family in the 2D lattice is restricted to those two gaps. Discrete solitons 
are associated with the T point of the first bancU^152l They have been observed in 
the context of optically-induced lattices by several group£ll 36 l 62 l Staggered solitons 
have also been observed in the case of defocusing nonlinearit}D^ and exist near the 
M-symmetry point of the first band. 

Gap solitons in self-focusing nonlinear media have been predicted to exist near 
the X-symmetry point at the top of the second bandSMl, having highly anisotropic 
structure and properties. These solitons are localized in one direction due to Bragg 
reflection and in the other due to total internal reflection. Their internal structure 
is fundamentally different from vortex solitons localized in the same gap but hav- 
ing symmetric profilesMES, or discrete solitons which can have elliptical profiles 
only due to their motion through the lattice 79 . To study the X-point gap solitons 
experimentall}^ we created a 2D square lattice and used a modulated extraor- 
dinarily polarized input beam [Fig. |9jb)] to match the profile of the Bloch wave 
associated with the X-symmetry point of the second band. At low laser power the 
beam diffracts and evolves into the corresponding Bloch wave [Fig.^c)]. Once the 
input laser power is increased, the output beam experiences a quasi-collapse to a 
reduced symmetry gap soliton as shown in Fig. [9jd,e). 

The important characteristic of this type of solitons is that they have highly 
anisotropic mobility properties. The solitons are mobile along their modulated x 
direction and highly immobile along the y direction. This mobility is illustrated 
when initial momentum is applied to the soliton. Our numerical simulations and 
experimental results are summarized in Fig. [9jf). An initial momentum along the 
y direction does not lead to any shift of the soliton, but only to a small modulation 
of its profile, similar to trapping of discrete solitons due to self-induced Peierls- 
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Fig. 9. Reduced symmetry gap solitons: In a two-dimensional square lattice (a) gap solitons 
can exist near the top of the second band (X-symmetry point), (b-e) Experimental excitation of 
reduced symmetry gap solitons: (b) input beam profile; (c) linear diffraction at the output; (d) 
soliton formation; (e) 3D representation, (f) Anisotropic soliton mobility when initial momentum 
is applied to the input beams. Left - numerical simulations; right - experimental results. 



Nabarro potential 85 l 86 l 87 l 88 l 89 l 90 l 91 l 9 ^lHl A momentum along x, on the other hand, 
leads to a shift of two lattice sites at the output. Our experimental observations 
well reproduced the numerical simulations [Fig. [9jf, right)]. The observed mobility 
properties of the reduced symmetry gap solitons relate to unique nonlinear transport 
of beams across the lattice, where the direction of beam propagation is determined 
by the localized state itself, and not by externally fabricated defects in the structure. 
The ability to robustly move along one particular direction of the lattice makes 
the reduced symmetry solitons good candidates for flexible soliton networks in 2D 
periodic structures. This offers new opportunities compared to the soliton networks 
suggested earlier for optical signal routing and switchin^^. 

Furthermore, the reduced symmetry gap solitons can be regarded as an analogue 
of a nonlinearity induced waveguide in periodic structures (for mobile solitons 
or as an optically induced high-Q cavity (for immobile solitons f^i. We believe that 
these ideas can be applied in other types of periodic structures such as photonic 
crystals and microstructured optical fibers. 



4. Fabricated periodic photonic structures 

Even though the optically-induced lattices offer great flexibility for dynamic mod- 
ification of the lattice parameters, experiments are complex and require extensive 
stability of the experimental setup in order to avoid fluctuation of the interfer- 
ence pattern of the lattice forming beams and averaging of the lattice potential. 
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Fig. 10. (a) Fabricated lithium niobate waveguide array, scheme of excitation, and intensity 
profile of the input beam, (b) Discrete diffraction at low laser power (10 nW): (top) 2D intensity 
profile, (middle) transverse intensity distribution and comparison with discrete model (crosses); 
(bottom) numerically calculated propagation inside the sample, (c) Formation of a staggered 
soliton at 1 mW laser power, soliton intensity profile and interferogram revealing the staggered 
phase structure, (bottom) Numerically calculated nonlinear propagation. 



A typical experimental setup covers roughly half an optical table, hence the direct 
applicability of the obtained results is limited. Therefore, it is important to trans- 
late the gained knowledge of nonlinear lattices into some more compact fabricated 
structures, offering closer connection to practical applications and devices. 



4.1. Discrete solitons in defocusing waveguide arrays 

To explore the propagation of light in fabricated periodic photonic structures, we 
utilized an array of optical waveguides fabricated by titanium indiffusion in a mono- 



crystal x-cut lithium niobate wafer [see Fig. 10 a, top)]. The array is 5cm long 
and consists of 100 closely spaced waveguides with a period of 19 or 10 fim. These 
fabricated structures offer great simplification of experiments. Furthermore, lithium 
niobate is known to exhibit strong photovoltaic nonlinearity of defocusing type^ at 
micro- Watt power levels for visible wavelengths. To test the linear propagation of 
light and demonstrate nonlinear localization in such arrays^, we coupled light into 
a single channel of the array as shown in Fig. [lO^a) and monitored the intensity 
profile at the output with increasing laser power. The input beam was circular with 
diameter of 2.7 am. At low laser powers (10 nW) the beam experiences discrete 
diffraction, where most of its energy is coupled away from the central waveguide 
[Fig.lIoFb)]. This behaviour is well reproduced by calculations using both a discrete 



model for weakly coupled lattice cites [crosses in Fig. 10 ^b, middle)], and a more 
general continuous model corresponding to our array (bottom). When the input 
laser power is increased light is coupled back into the central guide of the array, and 
at powers of about 1 mW all light is concentrated in the central guide of the array 
[Fig. [lO^c)]. An interferometric measurement of the output beam profile reveals 
that light in the central and the neighboring channels is out of phase, which is a 
direct indication that the observed localized state has a staggered phase structure 
and resembles a staggered soliton^. 
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The power dependent transition between discrete diffraction at low power where 
almost no light is transmitted inside the central channel, and the single channel 
localization proves to be beneficial for applications of laser mode locking^. In such 
application the spatial dynamics of light is translated to temporal dynamics and 
short pulse formation in the laser cavity. Furthermore, the waveguide technology 
is well compatible with fibers and applications of Kerr-lens mode-locking of fiber 
lasers leading to high contrast pulse generation seems feasible. The staggered phase 
structure of the localized state in the case of defocusing nonlinearity also appears 
to be more advantageous comparing to localization in structures with self-focusing 
nonlinear response^. We expect that the staggered profile can lead to improved 
contrast of the pulse generation. It is also important to note that the transition 
from diffraction to nonlinear localization has a well-defined threshold 97 , in contrast 
to the localization in homogeneous media. Such a threshold behaviour can be also 
applied to an optical diode device, replicating the transmission characteristics of 
usual electronic diodes. 



4.2. Nonlinear surface waves 

Additional functionality for control of light propagation in photonic lattices can be 
achieved when the optical beam interacts with a defect or an interface ( "surface" ) 
of the structure. For example, an optical beam injected at the edge waveguide 
of the array will be reflected from the interface formed between the periodic and 
homogeneous parts of the structure [Fig. pTj^a)]. This effective "repulsion" from 
the surface leads to modified discrete diffraction of the beamP^I and practically 
no light remains in the input edge waveguide. At high input powers, however, 
the surface can support various types of nonlinear surface waves 98 . In the case of 
focusing nonlinear response of the material, such nonlinear surface waves exist in 
the form of discrete surface solitons 98 . They have a propagation constant inside 
the total internal reflection gap of the periodic structure since the high power 
light induces a positive index defect at the edge waveguide of the array [Fig. 



11 



top)]. Such solitons were recently observed experimentally in AlGaAs waveguide 
arrays^. In the case of defocusing nonlinearity, surface solitons can also exist 98 . 
Their propagation constant resides inside the Bragg reflection gap and therefore 
they can be termed as surface gap solitons^°4 Due to the defocusing nonlinearity, 
the high intensity of the light results in a negative index defect at the edge waveguide 



[see Fig. 11 ^b, bottom)]. This defect leads to a less localized nonlinear surface wave 
which extends deeper inside the bulk medium. This implies that the surface gap 
solitons can be potential candidates for application of surface sensing and testing. 
Furthermore, unlike their discrete counterparts in self-focusing media, the surface 
gap solitons have a well-pronounced staggered phase structure, resulting in zero 
intensity in between the waveguides. 

In our experiments 101 we demonstrated the formation of surface gap solitons by 
injecting a narrow beam at the edge waveguide of a lithium niobate waveguide array. 
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Fig. 11. Surface waves: (a) Linear repulsion from the surface, (b) Theoretically predicted non- 
linear surface waves for focusing (top) and defocusing (bottom) nonlinearities. (c) Experimentally 
observed surface gap soliton in defocusing lithium niobate waveguide array. 



At low laser power (100 nW) the light is reflected by the surface [Fig. 11 'a, top)], 
while as the power is increased the beam is attracted to the surface and at a power 
of 0.5 mW, a surface gap soliton is formed in the array [Fig. 



11 



top)]. By use 

of interferometric measurements we also confirmed the staggered phase structure 
of the gap soliton. As seen in Fig. 



11 



'c, bottom) the interference fringes in each 
neighboring guide are shifted by half a period. 

Similar nonlinear surface waves have been demonstrated by other groups, 
in photorefractive iron-doped lithium niobate waveguide arrayd^l, and period- 
ically poled lithium niobate arrays with quadratic nonlinearit}A±£2J i n the latter 
case, the nonlinear response is dependent on the coupling between the fundamen- 
tal beam and the generated second harmonic, and thus the nonlinear response 
can be tuned from focusing to defocusing by temperature control, adjusting the 
phase mismatch of the harmonic generation process. Surface solitons may also 
be supported by optical nonlinearities with a nonlocal response^R More gen- 
eral types of nonlinear localized waves can exist at superlattice boundaries'^, in- 
terfaces between different lattices^ 06 ! 107 !, in two-dimensional geometries' 1Q8 I 1Q6 LL2£I, 
and at lattice dislocations' 11Q l m IHl The presence of interfaces and introduced 
or dynamically induced defects can also be used to perform tunable beam 



steering 



112|113|114|115|116|117|118| 



Additional possibilities are associated with multi- 



ple beam interactions and the formation of vectoiE5LLH2l Q r multi-garJ 5Q l 51 l 121 -HHI 
surface solitons. 



5. Concluding remarks 

In this review, we have described several fundamental effects observed in light prop- 
agation in nonlinear periodic structures. We would like to emphasis the ability of 
periodic structures to control the beam diffraction and refraction by tuning the 
structural parameters of the photonic lattice. We have also demonstrated how non- 
linear self-action of beams can be used to balance their diffraction through the 
formation of nonlinear localized states, or spatial optical solitons. These solitons 
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can have different properties and characteristics, depending of the type of nonlin- 
earity and their association with the specific Bloch waves of the periodic structure. 
Furthermore, we have demonstrated how the edges of the periodic structures can 
affect the nonlinear localization through the formation of nonlinear surface waves. 

In addition, we have also pointed towards possible applications of the observed 
effects. However, it is important to state that a wider range of applications is still 
to follow. We have shown that most effects are scalable and, therefore, they can be 
observed in structures with different periodicity. We were able to demonstrate the 
results in different periodic structures - from optically-induced lattices to fabricated 
waveguide arrays. The application of the presented concepts to smaller structures 
and their integration in photonic devices on an optical chip is an exciting direction 
to follow. We hope that with the current development of modern technologies and 
fabrication of photonic crystals a bright future of nonlinear periodic structures for 
light control on micro and nanoscales is soon to come. 
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